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Abstract 

Let X — {X(t),t 6 M. N } be a centered Gaussian random field with stationary in- 
crements and let T C 1* be a compact rectangle. Under X(-) G C 2 (K Ar ) and cer- 
tain additional regularity conditions, the mean Euler characteristic of the excursion set 
A u = {t G T : X(t) > u}, denoted by K{<p(A u )}, is derived. By applying the Rice 
method, it is shown that, as u — > oo, the excursion probability P{sup teT X(t) > u} 
can be approximated by E{</?(Au)} such that the error is exponentially smaller than 
E{(^(A„)}. This verifies the expected Euler characteristic heuristic for a large class of 
Gaussian random fields with stationary increments. 

Key Words: Gaussian random fields, stationary increments, excursion probability, excursion set, 
Euler characteristic, super-exponentially small. 

1 Introduction 

Let X = {X(t),t E T} be a real-valued Gaussian random field on probability space (O, J 7 , P), 
where T is the parameter set. The study for excursion probability P{sup tgr X(t) > u} 
is a classical but very important problem in probability theory and has many applications 
in statistics and related areas. Many authors have developed various methods for precise 
approximations of P{sup te <^ X(t) > u}. These include the double sum method [Piterbarg 
(1996a)], the tube method [Sun (1993)], the Euler characteristic method [Adler (2000), Taylor 
and Adler (2003), Taylor et al. (2005), Adler and Taylor (2007)] and the Rice method [Azais 
and Delmas (2002), Azai's and Wschebor (2005, 2008, 2009)]. 
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For a unit-variance smooth Gaussian random field X = {X(t),t £ T} parameterized on 
manifold T, Adler and Taylor (2007, Theorem 14.3.3) proved, under certain conditions on 
the regularity of X and topology of T, the following approximation: 



where <p(A u ) is the Euler characteristic of excursion set A u = {t £ T : X(t) > u} and 
a > is a constant which relates to the curvature of the boundary of T and the second order 
partial derivatives of X. This verifies the "Expected Euler Characteristic Heuristic" for unit- 
variance smooth Gaussian random fields. See also Piterbarg (1996a, 1996b), Takemura and 
Kurki (2002), Taylor and Adler (2003), Taylor et al. (2005) for similar results in special cases. 

The approximation (jl.ip is remarkable and very accurate, since E{(/9(^4 n )} is computable 
and the error is exponentially smaller than this principal term. It has been applied for 
P-value approximation in many statistical applications to brain imaging, cosmology, and 
environmental sciences. We refer to Adler and Taylor (2007) and its forthcoming companion 
Adler, Taylor and Worsley (2012) for further information. However the above requirement of 
"constant variance" on the Gaussian random fields is too restrictive for many applications and 
excludes some important Gaussian random fields such as those with stationary increments [see 
Section 2 below], or more generally, Gaussian random intrinsic functions [Matheron (1973), 
Stein (1999, 2012)]. If the constant variance condition is not satisfied, then the formulas for 
computing E{<^(j4 u )} [cf. Theorem 12.4.1 and Theorem 12.4.2 in Adler and Taylor (2007)] 
cannot be applied, and little had been known on whether the approximation (jl.ip still holds. 
The only exception is Azai's and Wschebor (2008, Theorem 5), where they proved (jl.ip for a 
smooth Gaussian random field whose maximum variance is attained in the interior of T. 

In this paper, let X = {X(t),t £ R^} be a centered real-valued Gaussian random field 
with stationary increments and let T C M> N be a rectangle. Our objectives are to compute 
E{9j(j4 u )} and to show that it can be applied to give an accurate approximation for the 
excursion probability P{sup tgr X(t) > u} and, in particular, (jl.ip holds for smooth Gaussian 
random fields with stationary increments. 

The paper is organized as follows. Firstly, in Section 2, we provide some preliminaries 
on Gaussian random fields with stationary increments and prove some basic lemmas. These 
lemmas are derived from the spectral representation of the random fields and will be useful 
for proving the main results in Section 3 and Section 4. 

In Section 3 we compute the mean Euler characteristic E{c^(^4 n )} by applying the Kac- 
Rice metatheorem in Adler and Taylor (2007, Theorem 11.2.1) [see also Adler and Taylor 
(2011, Theorem 4.1.1)]. The computation of E{(/?(j4 u )} involves the conditional expectation 
of the determinant of V 2 X(t) given X(t) and VX(t), which is more complicated for random 




as u — > oo 



(1.1) 
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fields with non-constant variance function. For Gaussian random fields with stationary in- 
crements, we are able to make use of the properties of VX and V 2 X to provide an explicit 
formula in Theorem 13.21 for E{(/j(^4 u )}, using only up to the second-order derivatives of the 
covariance function. 

Section 4 is the core part of this paper. Theorems 14.71 and 14.81 provide approximations 
to the excursion probability which are analogous to (jl.ip for Gaussian random fields with 
stationary increments. Since these random fields do not have constant variance, it is not 
clear if the original method for proving Theorem 14.3.3 in Adler and Taylor (2007) is still 
applicable. Instead, our argument is based on the Rice method in Azais and Delmas (2002) 
[see also Adler and Taylor (2007, p. 96-99)]. More specifically, we decompose the rectangle 
T into several faces of lower dimensions and then apply the idea of Piterbarg (1996b) and 
the Bonferroni inequality to derive upper and lower bounds for P{sup 4g *p X(t) > u} in terms 
of the the number of extended outward maxima [see (14. ip . (14. 2)) ] and the local maxima 
[see (|4.3|) ]. respectively. The main idea is to show that, in both cases, the upper bound 
makes the major contribution for estimating ¥{sup teT X(t) > u}, and the last two terms 
in the lower bounds in (|4.2p and (|4.3|) are super-exponentially small. Under an additional 
condition on the variogram of X, we apply (j4.3j) to obtain in Theorem 14. 71 an expansion for 
the excursion probability which is, in spirit, similar to the case of stationary Gaussian fields 
[cf. (14.0.3) in Adler and Taylor (2007)]. Theorem 14.81 establishes a general approximation to 
P{sup ter X(t) > u} in terms o{M{ip(A u )}, which verifies the "Expected Euler Characteristic 
Heuristic" for smooth Gaussian random fields with stationary increments. 

Section 5 provides further remarks on the main results and some examples where signif- 
icant simplifications can be made. In Example 15.31 and Example 15.41 we show that if the 
variance function of the random field attains its maximum at a unique point, then one can 
apply the Laplace method to derive a first-order approximation for the excursion probability 
explicitly. Finally, the Appendix contains proofs of some auxiliary lemmas. 

2 Gaussian Fields with Stationary Increments 
2.1 Spectral Representation 

Let X = {X(t),t E M. N } be a real- valued centered Gaussian random field with stationary 
increments. We assume that X has continuous covariance function C(t,s) = E{A(i)X(s)} 
and X(0) = 0. Then it is known [cf. Yaglom (1957)] that 

C(t, a) = [ {e l{t ' x) - l)(e" i(s ' A) - 1) F(d\) + (t, Os) (2.1) 
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where (x, y) is the ordinary inner product in M. N , is an N x N non-negative definite matrix 
and F is a non-negative symmetric measure on R^O} which satisfies 



/ 



1 + 1 A 



2 

F{dX) < oo. (2.2) 



Similarly to stationary random fields, the measure F and its density (if it exists) /(A) are 
called the spectral measure and spectral density of X, respectively. 

By (|2.ip we see that X has the following stochastic integral representation 

X(t) = f (e^> A > - l)W(dX) + (Y,t), (2.3) 
Jr n 

where Y is an A-dimensional Gaussian random vector and W is a complex- valued Gaussian 
random measure (independent of Y) with F as its control measure. It is known that many 
probabilistic, analytic and geometric properties of a Gaussian field with stationary increments 
can be described in terms of its spectral measure F and, on the other hand, various interesting 
Gaussian random fields can be constructed by choosing their spectral measures appropriately. 
See Xiao (2009), Xue and Xiao (2011) and the references therein for more information. 

For simplicity we assume that Y = 0. It follows from (|2.1|) that the variogram v of X is 
given by 

v(h) := R(X(t + h) - X{t)f = 2 / (1 - cos (h, X))F(dX). (2.4) 

Mean-square directional derivatives and sample path differentiability of Gaussian random 
fields have been well studied. See, for example, Adler (1981), Adler and Taylor (2007), Pot- 
thoff (2010), Xue and Xiao (2011). In particular, general sufficient conditions for a Gaussian 
random field to have a modification whose sample functions are in C k (W N ) are given by 
Adler and Taylor (2007). For a Gaussian random field X = {X(t),t G W N } with stationary 
increments, Xue and Xiao (2011) provided conditions for its sample path differentiability in 
terms of the spectral density function /(A). Similar argument can be applied to give the 
following spectral condition for the sample functions of X to be in C k (M. N ), whose proof is 
given in Cheng (2013) and is omitted here. 

Proposition 2.1 Let X = {X(t),t G W N } be a real-valued centered Gaussian random field 
with stationary increments and let ki (1 < i < N) be non-negative integers. If there is a 
constant £ > such that 



N 

/ n|A i | 2fc ^ e F(dA)<oo, (2.5) 

J{\\M\>i} f=i 
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then X has a modification X such that the partial derivative — f X ^\. is continuous on M N 
J F aq 1 -dt k N N 

almost surely, where k = Yli=i^i- Moreover, for any compact rectangle T C R N and any 
e' G (0,£ A 1), there exists a constant c\ such that 

/ d k x(t) d k x( s ) \ 2 „ UF , u 

E (— E L 4 i - L ^\ <cx\\t-s\\ E , Vt,sGT. (2.6) 

\dt k ^---dt k N N ds ki ■ ■ ■ ds k N N J ~ 

For simplicity we will not distinguish X from its modification X. As a consequence of 
Proposition [27T], we see that, if X = {X(t),t G R^} has a spectral density /(A) which satisfies 

fW = °( \\\\\ N + 2k+H ) as H A H ^ °°' ( 2 - 7 ) 

for some integer k > 1 and if 6 (0, 1), then the sample functions of X are in C k (WL N ) a.s. 
Further examples of anisotropic Gaussian random fields which may have different smoothness 
along different directions can be found in Xue and Xiao (2011). 

When X(-) G C 2 (R N ) almost surely, we write = X t (t) and §^ = De- 

note by VX(t) and V 2 X(t) the column vector (Xi(t), . . . ,XN(t)) T and the N x N matrix 
(Xij(t))ij=i t ,,, : N, respectively. It follows from (|2.1|) that for every t G 



:= / \ i \ j F(d\) = d £ ( -*' s) =E{Xi(t)Xj(t)}. (2.8) 

Let A = (Ajj)j^ = i jv, then (|2.8p shows that A = Cov(VX(t)) for all i. In particular, the 
distribution of VJ(t) is independent of t. Let 



Xij(t) := / AjAj cos (t, A) F(d\), A(t) := (A^i));,^!,...,^. 

JR^ 

Then we have 

Aii(*) - Aij- = f AA^cos <t, A) - 1) F(dX) = = E{X(t)X y (f)}, 

JgJV OtiOtj s=t 

or equivalents, A(t) - A = E{A(t)V 2 A(t)}. 

2.2 Hypotheses and Some Important Properties 

Let T = n£i[ a «' ^i] be a compact rectangle in M. N , where a% < hi for all 1 < i < N and ^ T 
[the case of G T will be discussed in Remark 15. 1| . In addition to the stationary increments, 
we will make use of the following conditions on X: 
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(HI). X(-) G C 2 (T) almost surely and its second derivatives satisfy the uniform mean-square 
Holder condition: there exist constants L and 77 G (0, 1] such that 

E(X^(i)-Xy(s)) 2 <L\\t-s\\ 2T >, Vi,sGT, i,j = l,...,N. (2.9) 
(H2). For every t £ T, the matrix A — A(t) is non-degenerate. 
(H3). For every pair (t, s) G T 2 with t ^ s, the Gaussian random vector 

(X(t),VX(t),X i:j {t), X(s),VX(s),X ij (s),l< i<j<N) 
is non-degenerate. 

(H3'). For every t G T, (I(t),VI(f),I i3 (t), 1 < i < j < N) is non-degenerate. 

Clearly, by Proposition 12. 1[ condition (HI) is satisfied if (|2.7p holds for k = 2. Also note 
that (H3) implies (H3'). We shall use conditions (HI), (H2) and (H3) to prove Theorems 
14.71 and 14.81 Condition (H3') will be used for computing E{i^(^4 n )} in Theorem 13.21 

The following lemma shows that for Gaussian fields with stationary increments, (H2) is 
equivalent to A — A(t) being positive definite. 

Lemma 2.2 For every t 7^ 0, A — A(t) is non-negative definite. Hence, under (H2), A — A(t) 
is positive definite. 

Proof Let t 7^ be fixed. For any (01, . . . ,a N ) G M^UO}, 

n r / N \ 2 

aiOjiXij - Xijit)) = j ( a { Xi J (1 - cos (t, A)) F(A). (2.10) 
i,j=i J RN ^ i=i ' 

Since ^h=\ a Ai) 2 (l — cos (f, A)) > for all A G R^, (|2.10p is always non-negative, which 
implies A — A(t) is non-negative definite. If (H2) is satisfied, then all the eigenvalues of 
A — A(t) are positive. This completes the proof. □ 

It follows from (|2.10p that, if the spectral measure F is carried by a set of positive Lebesgue 
measure [i.e., there is a set B C R^ with positive Lebesgue measure such that such that 
F(B) > 0], then (H2) holds. Hence, (H2) is in fact a very mild condition for smooth 
Gaussian fields with stationary increments. 

Lemma [2. 2 1 and the following two lemmas indicate some significant properties of Gaussian 
fields with stationary increments. They will play important roles in later sections. 

Lemma 2.3 Let t G M. N be fixed. Then for all i,j, k, the random variables Xi(t) and Xjf~(t) 
are independent. Moreover, M{Xij(t)Xki(t)} is symmetric in i, j, k, I. 
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Proof By (|2.ip . one can verify that for t, s G 



E{Xi(t)X jk (a)} = = / A * A A sin <* " s ' A > F ( dA )- 

Letting s = t we see that Xj(£) and Xj k (t) are independent. Similarly, we have 

EiXi, (t)X kl {s)} = jy% ,8 l = [ AiAjAfcA, cos (t - a, A) F(d\). 
dtidtjdskdsi J r n 

This implies the second conclusion. □ 

The following lemma is a directly consequence of Lemma 12.31 
Lemma 2.4 Let A = (o-ij)i<i,j<N be a symmetric matrix, then 

S t (i,j,k,l)=E{(AV 2 X(t)A) ij (AV 2 X(t)A) kl } 
is a symmetric function of i, j, k, I. 

3 The Mean Euler Characteristic 
3.1 Related Existing Results and Notation 

The rectangle T = n^Li[ a «> can be decomposed into several faces of lower dimensions. We 
use the same notation as in Adler and Taylor (2007, p. 134). A face J of dimension k, is defined 
by fixing a subset a(J) C {1, . . . , N} of size k and a subset e(J) = {£j,j ^ &(J)} C {0, l} N ~ k 
of size N — k, so that 

J = {t = (ti, . . . ,t N ) G T :aj < tj < bj if j G a( J), 

tj = (1 — £j)aj + Ejbj if j cr(J)}. 

Denote by d k T the collection of all /c-dimensional faces in T, then the interior of T is given by 

o 

T = djyT and the boundary of T is given by dT = U^Jq 1 U Je9fcT J. For J G c\.T, denote by 
VX|j(t) and V 2 X|j(i) the column vector (Xj^t), . . . , AQ fc (t))T ifegcr (j) and the fcxfc matrix 
(X mn (t)) mjne(T (j), respectively. 

If Al (•) G C 2 (K^) and it is a Morse function a.s. [cf. Definition 9.3.1 in Adler and Taylor 
(2007)], then according to Corollary 9.3.5 or page 211-212 in Adler and Taylor (2007), the 
Euler characteristic of the excursion set A u = {t G T : X(t) > u} is given by 

N k 

tp(A u )=j2 E (-i) fc E(- 1 )v<( J ) 

k=0 J€d k T i=0 
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(3.4) 



with 

Hi{J) := #{i G J : > u, VX\j(t) = 0, index(V 2 A| j(t)) = i, 
EjXjit) > for all j $ a(J)}, 

where e* = 2Ej — 1 and the index of a matrix is defined as the number of its negative 
eigenvalues. We also define 

Jli(J) := #{teJ: X(t) > u, VX {J (t) = 0, index(V 2 X| j(t)) = »}. (3.3) 

Let a 2 = Var(X(t)) and let = sup teT o\ be the maximum variance. For Gaussian 
fields with stationary increments, it follows from (j2.4|) that v(t) = a 2 . For t G J G <9fcT, 
where > 1, let 

Aj = {\j)i,jea{J)i A j(t) = (^'(*))j,ie(r(J)) 
2 = Var(X(t)|VX| J (t)), 7t 2 = Var(X(t)|VX(i)), 
{J 1 ,...,J N - k } = {l,...,N}\a(J), 
E{J) = {(tj,,. . .,tj N _ k ) G R N ~ k : ^e* > Q,j = Ji, Jiv-fc}- 

Then for all i G J, 

Aj = Cov(VX|j(*)), A J (t)-A J = E{X(t)V 2 X| J (t)}. (3.5) 

Note that Q\ > 7 2 for all t G T and 2 = 7i 2 if t G d^T. If {t} G d T, then VX|j(i) is 
not defined, in this case we set 6 2 as a 2 by convention. Let Cj(t) be the (1, j + 1) entry 
of (Cov(JT(t),VA'(t)))- 1 , i.e. 0,(4) = M lj+1 /detCov(X(t), VX(t)), where M 1J+1 is the 
cofactor of the (l,j + 1) entry, K{X(t)Xj(t)}, in the covariance matrix Cov(X(t), VX(t)). 

Denote by H k (x) the Hermite polynomial of order fc, i.e., H k (x) = (—l) k e x2 ^ 2 -^(e~ x2 ^ 2 ). 
Then the following identity holds [cf. Adler and Taylor (2007, p. 289)]: 

/•oo 

/ H k (x)e- x2 / 2 dx = H k ^(u)e- u2 / 2 , (3.6) 

J u 

where u > and fc > 1. For a matrix A, we will use \A\ or detA to denote its determinant. 
Let M+ = [0,oo), M_ = (-oo,0] and *(«) = (2vr)- 1 /2 f™ e ~ x2 / 2 dx. 

3.2 Computing the Mean Euler Characteristic 

The following lemma is an analogue of Lemma 11.7.1 in Adler and Taylor (2007). It provides 
a key step for computing the mean Euler characteristic in Theorem 13.21 meanwhile, it has 
close connection with Theorem 14.71 
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Lemma 3.1 Let X = {X(t),t G M. N } be a centered Gaussian random field with stationary 
increments satisfying (HI), (H2) and (H3'). Then for each J € dj-T with k > 1, 

-{ h-^} - (2 J-^r- 1 0.7, 

Proof Let T>% be the collection of all k X k matrices with index i. Recall the definition 
of Jli(J) in (|3.3p . thanks to (HI) and (H3'), we can apply the Kac-Rice metatheorem [cf. 
Theorem 11.2.1 or Corollary 11.2.2 in Adler and Taylor (2007)] to get that the left hand side 
of (|3.7p becomes 

/ Pvx |jW (0)(it^(-l) i E{|detV 2 X| J (t)|l {v2X|j(t)eI , i} l {xw > M} |VX| J (t) = 0}. (3.8) 
Jj i=0 

Note that on the event Pj, the matrix V 2 A| j(i) has i negative eigenvalues, which implies 
(-l) i |detV 2 X| J (t)| = detV 2 X|j(t). Also, uf =0 {V 2 X|j(t) £Dj = fi a.s., hence $£E§ equals 

Pvx |jW (0)^E{detV 2 A| J (t)l { x W > u} |VX| J (t) = 0} 

f e -x 2 /(29?) roo ( 3 - 9 ) 

= Jj (2n)(^\Aj\^e t dt I dxE{deW% J (t)\X(t)=x,VX lJ (t) = 0}. 

Now we turn to computing E{detV 2 A| j(t)\ X(t) = x,VX\j(t) = 0}. By Lemma 
under (H2), A — A(t) and hence Aj — Aj(i) are positive definite for every t £ J. Thus there 
exists a k x k positive definite matrix Qt such that 

Q t (Aj - Aj(t))Q t = I k , (3.10) 

where is the k x k identity matrix. By ()3.5|) . 



E{A(t)(Q 4 V 2 A| J (t)Q t ) 4 ,} = -(Q t (Aj - Aj(t))Q t ) fj = -<%, 

where <5jj is the Kronecker delta function. One can write 

E{det(Q t V 2 A| j(t)Q t )| A(t) = z, VA, j(t) = 0} = E{detA(t, x)}, (3.11) 

where A(t,x) = (Ajj(t, x))j JgCT (j) with all elements Aij(t,x) being Gaussian variables. To 
study A(t,x), we only need to find its mean and covariance. Note that VA(t) and V 2 A(i) 
are independent by Lemma 12.31 then we apply Lemma 16.11 to obtain 

E{Ay(t,x)} =E{(Q t V 2 X l j(t)Q t ) lJ \X(t) = x,VX ]J (t) = 0} 

= (K{X(t)(Q t V 2 X l j{t)Q t ) l3 ], 0, . . . , 0)(Cov(X(t), VA, j(t)))~\^ 0, . . . , 0) T (3 . 12) 

= (Sij, 0, . . . , 0)(Cov(A(t), VXuCt)))- 1 ^, 0, . . . , Of = --J<%, 
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where the last equality comes from the fact that the (1,1) entry of (Cav(X(t), VX| j(t))) -1 
is detCov(VX|j(i))/detCov(X(t), VX|j(t)) = l/0f. For the covariance, applying Lemma IBTTl 
again gives 

E{(A ii (t, x) - E{A y (t, x)})(A w (t, x) - E{A fci (t, x)})} 

= E{(Q t V 2 X| j(t)Q t )y (Q t V 2 X| j(i)Qt) fc/ } - (E{X(t)(Q t V 2 *|.,(t)Q t )ii}> 0, . . . , 0) 

• (Cov(X(t), VX| J (t)))- 1 (E{X(t)(Q i V 2 X| J (t)Q 4 ) M }, 0, . . . , 0) T 
= j, k, I) - (-<%, 0, . . . , 0)(Cov(X(i), VX^m^i-Sku 0, . . . , 0) T 

&ij&kl 

= S t {i,j,k,l) , 

where St is a symmetric function of i, j, k, I by applying Lemma 12.41 with ^4 replaced by Qt- 
Therefore (|3.1ip becomes 
1 



E<{ ^det(^Q t (V 2 Jf| J (t))Q t ) X(t) = x,VX {J (t) = } = ^E<j det( A(i) - f 4 



^{det(z 



x 



where A(i) = (A^(f))j t j^ a (J) an d an Ajj(i) are Gaussian variables satisfying 

E{Ay(i)} = 0, E{Aij{t)A kl {t)} = efS t {iJ,k,l) - 5 i:i 6 k i. 
By Corollary 11.6.3 in Adler and Taylor (2007), (|3TTT) is equal to (-l) k e^ h H k (x/e t ), hence 
E{detV 2 X|j(t)|A(t) = x,VX\j{t) = 0} 

= E{<let{Q^Q t V 2 X\j{t)Q t Q^)\X{t) = x, VX|j(i) = 0} 
= \Aj - Aj(t)|E{det(Q t V 2 X|j(t)Q t )|X(t) = x, VX |i7 (t) = 0} 
(-1)' 



IAj-Aj^IF, 



't 



Plugging this into (|3.9p and applying (j3.6j) . we obtain the desired result. □ 

Theorem 3.2 Let X = {X(t),t G R^} &e a centered Gaussian random field with stationary 
increments such that (HI), (H2) and (H3') are fulfilled. Then 

N 1 

%^}= E nx(t)>u,vx(t)eE({t})) + i: E faww 

{Oea T fc=i Jed k T v j 1 J 1 

f , f°° , f f r , |Aj-Aj(t)| 

x #fc + 7*Cji (*)l/Ji + • • • + 7tCj*_* (t)yj N - k ) 

x Px(t),x Jl (t),...,x Jjv _ fc ( i )(z,2/Ji, • • • ,2/J^_jVA|j(t) = 0). 
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Proof According to Corollary 11.3.2 in Adler and Taylor (2007), (HI) and (H3') imply 

that X is a Morse function a.s. It follows from (|3.1|) that 

N f k , 

EMA U )} = Y,H ("W D" 1 ^^) • (3-14) 
k=0 J£d k T i=Q ' 

If J € 9 T, say J = {t}, it turns out E{fi (J)} = P(X(t) > u,VX(t) G £({£})). If J G d k T 
with fc > 1, we apply the Kac-Rice metatheorem to obtain that the expectation on the right 
hand side of (|3.14p becomes 

^vx |jW (0)dt^(-l) i IE{|detV 2 A| 



x l{x( t )> u} |VX| J (t)=0} 

(2^/^,1^ 7j -J u dX J-J E{J) d ^- dy ^ (3 ' 15) 
x E{detV 2 X|j(i)|X(t) = x, Xj.it) = VJl ,.. .,Xj N _ h (t) = yj N _ k ,VX\j(t) = 0} 
x Px{t),x Jl {t),...,Xj lf _ k ( L t)(x,yj 1 ,. ■ ■ ,yj N _ k \VX VJ {t) = 0). 

For fixed t, let Qt be the positive definite matrix in (|3.10|) . Then, similarly to the proof in 
Lemma we can write 

E{det(Q t V 2 X l j(t)Q t )\X(t)=x,X Jl (t) = yj 1 ,...,Xj N _ k =yj N ^,VX l j(t) = 0} 
as E{detA(i, x)}, where A(t, x) is a matrix consisting of Gaussian entries Ay(t, x) with mean 

E{(Q t V 2 X|.,(t)Qi)yl*(*) = x,X Jl (t) = yjl , . . . , A Jjv _ fc = y Jjv _ fc , VA|j(t) = 0} 
= (-<%, 0, . . . , 0)(Cov(A(i), X, 7l (i), . . . , Xj N _ k (t), VA, At)))- 1 

■(x, yjl ,...,yj N _ k ,0,...,0) T ( 3 - 16 ) 

= --^(x + llC Jx {t) yjl +■■■ + i 2 t Cj N _ k (t)yj N _ k ), 
and covariance 

E{ (Aij (t,x)- E{ A^ (i, x)})(K H (t,x)- E(A kl (t,x)})} = S t (i,j,k,l) - . 

Following the same procedure in the proof of Lemma 13. 1| we obtain that the last conditional 
expectation in (|3. 15[) is equal to 

(-l) k \Aj- Aj(t)\ H ^ + 7tCji(t)yji + . . . + ltC j N „ k (t)yj N _ k ) . (3.17) 
Plug this into f)3. 15H and (|3.14p . yielding the desired result. □ 
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Remark 3.3 Usually, for nonstationary (including constant-variance) Gaussian field X on 
R , its mean Euler characteristic involves at least the third-order derivatives of the covari- 
ance function. For Gaussian random fields with stationary increments, as shown in Lemma 
12.31 "&{Xij(t)Xk(t)} = and E{Xy (i)-Xjy(i)} is symmetric in i, j, k, I, so the mean Euler 
characteristic becomes relatively simpler, contains only up to the second-order derivatives 
of the covariance function. In various practical applications, (|3.13[) could be simplified with 
only an exponentially smaller difference, see the discussions in Section 5. 

4 Excursion Probability 

4.1 Preliminaries 

As in Section 3.1, we decompose T into several faces as T = \J^ =Q dkT = (JfcLo Uje<9 fc T J- 
For each J € <9/%T, define the number of extended outward maxima above level u as 

Mf (J) := #{t€ J: X(t) > u,VX\j(t) = 0, index(V 2 X| j(t)) = k, 

e*Xj(t) > for all j $ a(J)}. 

In fact, M^{J) is the same as /Xfe(J) defined in (|3.2|) with i = k. We will make use of the 
following lemma, whose proof is given in the Appendix. 

Lemma 4.1 Let X = {X(t),t € M. N } be a Gaussian random field satisfying (HI) and (H3') ; 
then for any u > 0, 

TV 

{ supX(t) > u} = |J |J {M* (J) > 1} a.s. 

teT k=0 J&d k T 

It follows from Lemma |4. II that 

TV TV 

p{su P a^) >^}<e E nM*(j) > i} < E E E i M u(J)}- (4.i) 

teT fc=o JGO fc T fc=o Jea fe T 

On the other hand, by the Bonferroni inequality, 

TV 

p{supX(t) >^}>E E F i M u(J) > 1} - E F i M u(J) > l,Mf (J') > 1}. 

t£T k=0 j£d k T J+J' 

Let Pi = F{M*(J) = i}, then P{Mf(J) > 1} = TZiPi and E{Mf(J)} = ££i*Fi> it 
follows that 

oo 

E{M* (J)} - P{Mf (J) > 1} = £(i - l) ft 

i=2 

< E = ( J )« ( J ) - 1)}- 
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Together with the obvious bound P{M^(J) > 1, M®(J') > 1} < E{M^(J)M^(J')}, we 
obtain the following lower bound for the excursion probability, 



»{supX(t) > u} > J2 E_ ( E {<( J )} - /HMf (J)(Mf (J) - 1)} 

(4-2) 

- £ E{Mf (J)Mf (J')}. 

Define the number of local maxima above level u as 

M U (J) := #{t£J: > u,VX |i7 (t) = 0, index(V 2 X| = A;}, 

then obviously M U (J) > Af^(J), and M„(J) is the same as fJ-k(J) defined in (I3.3P with i = k. 
It follows similarly that 



Yl E i M u(J)} > P{ supX(t) > u} 

«w - (43) 

> J] ^ ( E { M «(^)} - ]fWu{J)(M u {J) - 1)}] - Y E{M U (J)M U {J')}. 

We will use (|4.ip and (|4.2p to estimate the excursion probability for the general case, 
see Theorem 14.81 Inequalities in (|4.3p provide another method to approximate the excursion 
probability in some special cases, see Theorem l4,7[ The advantage of (|4.3h is that the principal 
term induced by Ylk=Q Sjea t E {-^"(^)} * s mucn easier to compute compared with the one 
induced by £f=o E W E{M* (■>)}• 

4.2 Estimating the Moments: Major Terms and Error Terms 

The following two lemmas provide the estimations for the principal terms in approximating 
the excursion probability. 

Lemma 4.2 Let X be a Gaussian field as in Theorem \3.2l Then for each J £ d^T with 
k > 1, there exists some constant a > such that 

■KM.W) - / ^^^.(^e-VW) dt(1 + (4 . 4) 
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Proof Following the notation in the proof of Lemma 13. 11 we obtain similarly that 
E{M U (J)} = ^pvx |J (t)(0)dtE{|detV 2 X| J (t)|l {V 2 X|j(i)e23fe} l {x(t) > u} |VX| J (t) = 0} 

r poo / j\& — a; 2 /(2f? 2 ) 

= J/ 1 L dX (2^+^|A J |V^, E{detV2X ' j(t)1 ^ v2x i^W gP ^ |X(t) = X ' VX \ j{t) = ° } - 

(4.5) 

Recall V 2 X\j{t) = Q^QtV 2 '■X\ J {t)Q t Q~[ 1 and we can write (137121) as 



E{Q t V 2 X {J (t)Q t \X(t) = x,VX {J (t) = 0} = --j/fc. 



.7.- 

Make change of variables 

V(t) = Q t V 2 X { j(t)Qt + -J/ fc , 

where V(t) = (Vjj(i))i<jj<fc. Then (V(t)\X(t) = x, VXij(i) = 0) is a Gaussian matrix whose 
mean is and covariance is the same as that of (QtV 2 X\ j(t)Qt\X(t) = i,Vlu(t) = 0). 
Denote the density of Gaussian vectors ({V%j(t))i<i<j<.k\X(t) = x, VX\j(t) = 0) by ht(v), 
v = («y)i<i<j<* e M fc ( fc+1 )/ 2 , then 



E{det(Q t V 2 X |J (i)Q t )l {V 2 X|j(t)eI?fc} |X(t) = x, VX|j(t) = 0} 

= E{det(Q t V 2 X, J (t)Q t )l {QtV2X|jWQteI?fc} |X(t) = x, VX {J (t) = 0} 

= / det ( (vij) - -^I k ) h t (v) dv, 

t 

where (u^) is the abbreviation of matrix (%)i<ij<fc. Since {# 2 : i G T} is bounded, there 
exists a constant c > such that 



K-)-|,4eP fc , v||( % -)lh= f E 4) 



1/2 

' X 

< -. 

c 



Thus we can write (14.611 as 



/ detf (vij) - -^4 )ht(v)dv - / detf (%•) - ^J fc ) ^(u) (fo 



(4.7) 



E{det(Q t V 2 X|j(t)Q t )|X(i) = x,VX\j{t) = 0} + Z(t, 



where x) is the second integral in the first line of (|4.7p and it satisfies 



|Z(t,x)| < 



:»«)ii>! 



X 
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ht{v)dv. 



Denote by G(t) the covariance matrix of ((Vij(t))i<i<j<k\X(t) = x,VX\j(t) = 0), then by 
Lemma 16.21 in the Appendix, the eigenvalues of G(t) are bounded from above and below by 
positive and finite constants for all t £ T, hence so are those of (G(i))" 1 for all t 6 T. It 
follows that there exists some constant a' > such that ht(v) = o(e _a '"^ iJ ^' 2 ) and hence 
= o(e~ ax ) for some constant a > uniformly for all t € T. Combine this with 
(|4.5p . (|4.6p . (|4.7p and the proof of Lemma [3TTI yielding the result. □ 

Lemma 4.3 Let X be a Gaussian field as in Theorem \S.°A Then for each J G d^T with 
k > 1, £/iere exists some constant a > suc/i i/iaf 



e{<(j)} = I dt /> / • • • L Avh ■ ■ ■ dyj 



x ' Aj f j(t)l ^ fc f - + ltG h (t) yjl + ■■■ + lt C jN _ k (t)y jN _ . \ 



7? V7t 



xPx(t),x Jl (t),...,x Jjv _ jfe (t)(^2/Ji,---,yj JV _JVX| J (f) = 0)(l + o(e Q " )). 

Proof Similarly to the proof in Theorem 13.21 we see that E{Mf (J)} is equal to 
(-l)*|Aj-Aj(t)| 



/ ( (zlPlA^i/i*^ /f dX / • • • / £(J) dy - h • • • ^-^{det(Q t V 2 X, j( t)Q t ) 

x 1 {Qt^ 2 x lJ (t)Q t ev k }\X(t) =x,X Jl (t) =y Jl ,...,Xj N _ k (t) = yj N _ k , VX\j{t) = 0} 

x Px(t),x. h (t),...,Xj N _ k (t) (x, 2/Ji , • • • , I VX| j(t) = 0) 
(-l)*|Aj-Aj(t)| 



where is the positive definite matrix in (|3.10p . Then using similar argument in the proof 
of Lemma 14.21 to estimate K(t,x), we obtain the desired result. □ 

Denote by S n C M n+1 the n-dimensional unit sphere. We call a function h(u) super- 
exponentially small [when compared with P(sup tgr X(t) > u)], if there exists a constant 
a > such that h(u) = o(e~ au2 ~ u ' ! /( 2o f)) as u — > oo. 

The following lemma is Lemma 4 in Piterbarg (1996b). It shows that the factorial mo- 
ments are usually super-exponentially small. 

Lemma 4.4 Let {X(t) : t £ M, N } be a centered Gaussian field satisfying (HI) and (H3). 
Then for any e > 0, there exists E\ > such that for any J £ dkT and u large enouth, 

E{M U (J){M U (J) - 1)} < e -« 2 /(2/3 2 7 + £ ) + e -u 2 /(2aj-e 1 )^ 

where f3j = sup tG j sup eG§ fe-i Vax(X(t)\VX\ j(t), V 2 X\ j(t)e) and ex 2 = sup tGj Var(X(t)). 
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Corollary 4.5 Let X = {X(i),t £ R^} be a centered Gaussian random field with stationary 
increments satisfying (HI), (H2) and (H3). Then for all J £ <9 fc T, E{M U (J)(M„(J) - 1)} 
and E{M,f (J)(M^ (J) — 1)} are super- exponentially small. 

Proof Since (J) < M W (J), we only need to show that E{M U (J)(M U (J) - 1)} is super- 
exponentially small. If k = 0, then M U ( J) is either or 1 and hence E{M M ( J)(M U ( J) — 1)} = 
0. If k > 1, then, thanks to Lemma 14.41 it suffices to show that 0j is strictly less than a\. 
Clearly, Var(A(t)|VA| j(t), V 2 X\ j(t)e) < a\. Applying Lemma O yields that 

Yw(X{t)\VX\j{t)y 2 X\j(t)e) = 4 => E{X(t)(V 2 X|j(t)e)} = 0. 

Note that the right hand side above is equivalent to (Aj(i) — Aj)e = 0. By (H2), Aj(t) — Aj 
is negative definite, which implies (Aj(t) — Aj)e / for all e £ S fc_1 , so that 

sup Vai(X(t)\V\jX(t),VfjX(t)e) <o\. 

Therefore /3 2 < cr 2 -, by continuity. □ 

The following lemma shows that the cross terms in (|4.2p and (|4. 3[) are super-exponentially 
small if the two faces are not adjacent. For the case when the faces are adjacent, the proof 
is more technical, see the proofs in Theorems 14.71 and 14.81 

Lemma 4.6 Let X = {X(t),t £ W N } be a centered Gaussian random field with stationary 
increments satisfying (HI) and (H3). Let J and J' be two faces ofT such that their distance 
is positive, i.e., inf 4e j jSg j/ \\s — t\\ > So for some So > 0, then E{M U ( J)M U { J')} is super- 
exponentially small. 

Proof We first consider the case when dim(J) = k > 1 and dim(J') = k' > 1. By the 
Kac-Rice metatheorem for higher moments [the proof is the same as that of Theorem 11.5.1 
in Adler and Taylor (2007)], 

E{M u (J)M„(JO} = ^di^a' S E{|detV 2 A| J (t)||detV 2 A| J K S )|l {xw > UiX{s) > u} 

x Mv*X\ J (t)eD k ,v'X\j,(8)eD k ,}\ x ( t ) = x > x ( s ) = V> vx \j( t ) = °> 
v ^|j'( s ) = ^}Px(t),x(s),vx l j(t),vx lJ ,( s ){x,y,0,0) 

r f f'OO poo 

< dt ds dx dyE{|detV 2 A|j(t)||detV 2 A|j,(s)| 

J J J J' Ju J u 

\X(t) = x,X{s) = y,VX\j(t) = 0, VX|j,(s) = 0}p X (t),X(s)(x,y) 
x Pvx lJ (t),vx, JI (s){0,0\X{t) =x,X(s) = y). 



(4.9) 
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Note that the following two inequalities hold: for constants a% and bj, 

k k' / k k' 



nwnfeis^(Ew +t '+EfeH' 

8=1 j=l V i=l j=l 7 

and for any Gaussian variable £ and positive integer I, 

E\C\ l <E(|E£| + |£-E£|/ < 2'(|E£|' + E|£-E£|') < 2 l (\E£\ l + C z (Var(e)) //2 ), 

where the constant C; depends only on /. Combining these two inequalities with Lemma l6.1| 
we get that there exist some positive constants C\ and N\ such that for large x and y, 

sup E{|detV 2 X| J (t)||detV 2 X| J ,(s)||X(t) = x,X{s) = y, 
ieJ,seJ' (4-10) 

VX| 7 (t) = 0, VX|j,( S ) = 0} < d^y**. 

Also, there exists a positive constant C*2 such that 

sup PvXij^vx^^OlX^) = x,X(s) = y) 
teJ,seJ' ' 

< sup (2vr)-( fc+fc ')/ 2 [detCov(VX J (t), VXj/(s)|X(t) = x,X(s) = y)]" 1 / 2 < C 2 . 
tej,s£j' 

Let jo(^o) = sup|| s _ t || >( s M^^Ifll w hich is strictly less than 1 due to (H3), then Ve > 0, 
there exists a positive constant C3 such that for all t E J, s S J' and it large enough, 

00 roo 



x Nl y Nl p X (t),x(s)(x,y)dxdy = E{[X(t)X(s)] N ' l {X {t)>u,x( s )>u}} 

( u 2 \ (412) 

< E{[X(t) +X(s)] 2N n {x{t)+x{s) y 2u} } < C 3 exp (eu 2 - j . 

Combine gSJ with P~TUD . (fl~TTD and (fl~T2"D . yielding that E{Mf (JjAff (J ; )} is super- 
exponentially small. 

When only one of the faces, say J, is a singleton, then let J = {t } and we have 



00 ^00 



E{M U (J)M U (J')} < / ds / dx / dy Px{to)>x{s)jVX {s) {x,y,0) 

J J' Ju Ju 1 (4.13) 

xK{\detV 2 X l j l (s)\\X(t ) = x,X(s) = y,VX lJ ,(s) = Q}. 

Following the previous discussion yields that E{M U (J)M U (J')} is super-exponentially small. 

Finally, if both J and J' are singletons, then E{M U ( J)M U { J')} becomes the joint proba- 
bility of two Gaussian variables exceeding level u and hence is trivial. □ 
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4.3 Main Results and Their Proofs 

Theorem 4.7 Let X = {X(t) : t £ M. N } be a centered Gaussian random field with stationary 
increments such that (HI), (H2) and (H3) are fulfilled. Suppose that for any face J, 

{t £ J : v{t) = c4, Uj(t) = for some j <£ a (J)} = 0. (4.14) 

Then there exists some constant a > such that 

N 



»\supX(t) >u\ = V V E{M U (J)} +o[e 



■au 2 -u 2 /(2aj,)\ 



teT k=0 J£d k T 



N i 

{t hr [J kjkT^ {k+1)/2 \^ 1/2 



J °t 



0$ 



Proof Since the second equality in ()4.15|) follows from Lemma 14.21 directly, we only need 
to prove the first one. By (|4.3p and Corollary 14.51 it suffices to show that the last term in 
(|4.3p is super-exponentially small. Thanks to Lemma 14.61 we only need to consider the case 
when the distance of J and J' is 0, or I := J D J' 7^ 0. Without loss of generality, assume 

a (J) ={l,...,m,m+l,...,k}, a(J') = {1, . . . , m, k + 1, . . . , k + k' - m}, (4.16) 

where 0<m<k<k'<N and k' > 1. If k = 0, we consider <r(J) = by convention. Under 
such assumption, J S <9fcT, J' G Sfc'T and dim(I) = m. 

Case 1: = 0, i.e. J is a singleton, say J = {to}. If ^(io) < <?r> then by (|4.13p . it is 
trivial to show that E{M U ( J)M U (J')} is super-exponentially small. Now we consider the case 
v(to) = a\. Due to (|4.14p . K{X(to)Xi(to)} 7^ an d hence by continuity, there exists 5 > 
such that E{X(s)X 1 (s)} / for all ||s-i || < 6. It follows from (|4TT3|) that E{M U (J)M U (J')} 
is bounded from above by 

/ ds dx dy E{|detV 2 X\ j> (s)\\X (t ) = x, X(s) = y,VX l j, (s) = 0} 

J s£j':\\s— to\\>S Ju Ju 

x Px(t () ),x(s),vx ]jl (s)(x,y,0) 

f' f'OO 

+ ds dyK{\detV 2 X lJ ,(s)\\X(s)=y,VX l Xs) = 0}p x{s) y X . Jl{s) (y,0) 

J s£j':\\s-to\\<6 Ju 

:= h + h- 
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Following the proof of Lemma 14.61 yields that I\ is super-exponentially small. We apply 
Lemma |6. II to obtain that there exists £q > such that 

sup Var(X(s)|VX|j,(s)) < sup Var(X(s)|X x (s)) < a\ - e Q . 

seJ':||s-to||<<5 seJ':||s-to||<<5 

Then I2 and hence E{M U ( J)M U (J')} are super-exponentially small. 

Case 2: k > 1. For all t G I with v(t) = of, by assumption flUH) , E{X(t)Xi(t)} ^ 0, 
V i = m + 1, . . . , k + k' — m. Note that I is a compact set, by Lemma 16.11 and the uniform 
continuity of conditional variance, there exist £i,5i > such that 

sup Vax(X(t)\X m+1 (t),...,X k (t),X k+1 (s),...,X k+k ,„ m (s)) < o\ - e u (4.17) 

t£B,s£B> 

where B = {t G J : dist(i, J) < <5i} and B' = {s £ / : dist(s,I) < <5i}. It follows from P~9"|) 
that E{M U (J)M U (J')} is bounded by 

/* /* /*oo roc 

/ / / dx dyp X (t),x(s),vx u (t),vx n/ (s){x,y,0,0) 

J J{JxJ')\(BxB>) Ju Ju 1 

x E{|detV 2 X|j(t)||detV 2 X|j,(s)||X(t) = x,X(s) = y,VX\j(t) = 0,VX| j,(s) = 0} 

+ // dtds / (tep A:(t) (x|V^|j(t) = 0,VX|j,(s) = 0Kx J (i),vx IJ ,( s )(0,0) 

x E{|detV 2 X|j(t)||detV 2 X|j,(s)||X(t) = a?, VX |i7 (t) = 0,VX\j>(s) = 0} 
:= / 3 + h- 
Note that 

(J x J')\(B x S') = (( J\B) x 5') |J (B x ( J\B)) |J (( J\B) x ( J\B)) . (4.18) 

Since each product set on the right hand side of (|4.18p consists of two sets with positive 
distance, following the proof of Lemma I4T61 yields that I3 is super-exponentially small. 
For I4, taking into account (|4.17p . one has 

sup yax(X(t)\VX\j(t),VX\ JI (8))<<4-e 1 . (4.19) 

To estimate 

Pvx^vX^s^O) = (2vr)-( fc+fe ')/ 2 (detCov(VX| J (t),VX| J ,( S )))- 1/2 , (4.20) 

we write the determinant on the right hand side of (|4.20p as 

detCov(X m+1 (t), . . .,X k (t),X k+1 (s), . . . (s)|Xi(t), . . . ,X m (t),Xi(s), . . .,X m (s)) 

x detCov(Xi(i), . . . ,X m (t),X 1 (s), . . . ,X m (s)), 

(4.21) 
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where the first determinant in (|4.2ip is bounded away from zero due to (H3). By (HI), as 
shown in Piterbarg (1996b), applying Taylor's formula, we can write 

VI(s) = VX(t) + V 2 X(t)(s -t) T + \\s- t\\ 1+ W ttS , (4.22) 

where Y tjS = (Y^ s , . . . , Y^) T is a Gaussian vector field with bounded variance uniformly for 
all t G J, s G J'. Hence as \\s — t\\ — > 0, the second determinant in (|4.21|) becomes 

detCov(Xi(t), . . . ,X m (t), Jfi(t) + (VXi(t), s -t) + \\s- t\\ 1+r >Y t ] s , . . . , 
X m (t) + (VX m (t),s -t) + \\s - t\\ 1+ ^Y t m s ) 
= detCov(Xi(t), . . . , X m (t), (VXx(t), s - t) + ||s - if^Y 1 ,, . . . , (4.23) 

(vx m (t), s -t) + || s -t|| 1+ "y™) 

= ||s - i|| 2 ™detCov(Xi(t), . . . ,X m (t), (VX 1 (t),e t>s ), (VX m (i), e M ))(l + o(l)), 

where e< )S = (s — t) T /||s — 1\\ and due to (H3), the last determinant in (|4.23p is bounded away 
from zero uniformly for all t G J and s £ J'. It then follows from (|4.2ip and f|4.23|) that 

detCov(VX|j(t), VX|j,(s)) > C x \\a - t\\ 2m (4.24) 
for some constant C\ > 0. Similarly to (|4.10p . there exist constants C 2 ,Ni > such that 
sup E{|detV 2 X|j(t)||detV 2 X|j,(s)||X(t) = x,VX\j(t) = 0, VX\j,(s) = 0} 
<C 2 (l + x Nl ). 

Combining (j4TTU|) with (]4T20|) . flOj]) and (|4T25]> . and noting that m < k! implies l/||s - t\\ m 
is integrable on J x J', we conclude that I4 and hence M{M U (J)M U (J')} are finite and super- 
exponentially small. □ 

Theorem 4.8 Let X = {X(t) : t G M. N } be a centered Gaussian random field with stationary 
increments such that (HI), (H2) and (H3) are fulfilled. Then there exists some constant 
a > such that 

N 

P{ supX(t) > u\ = V y E{Mf (J)} + ( e -^ 2 -«7(2<4)) 
lieT J ' (4-26) 

= E{(^(^)} + (e-^ 2 - u2 /( 2 4)), 

where E,{ip(A u )} is formulated in Theorem \3.S\ 
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It is worth mentioning here that the main idea for the proof of Theorem 14.81 comes from 
Azai's and Delmas (2002) [especially Theorem 4]. Before showing the proof, we list the 
following two lemmas whose proofs are given in the Appendix. 

Lemma 4.9 Under (H2), there exists a constant Qo > such that 

(e, (A - A(i))e) > a , V t € T, e € S^" 1 . 
Lemma 4.10 Let {£i(t) : t £ T%} and {^(t) '■ t G T2} be two Gaussian random fields. Let 
crf(i) = Var(£j(f)), o=j = super^i), 0^ = inf cr^t), 



teT, 



anc? assume < a, < <fj < 00, where i = 1,2. If < p < p < 1, then for any N±, N2 > 0, 
there exists some a > smc/i i/iai as u ^ 00, 

snp E{(1 + + k 2 ( S )r 2 )l {?1 ( t )>^ 2 W<0}} = o(e— 2 -« 2 /^)). 

t£T 1 ,seT 2 

Similarly, if — l<p<p<0, t/ien 

snp E{(1 + |&(t)|^ + |6( S )r 2 )l {?1 W>^ 2 (s)>0}} = o(e— 2 -" 2 /^)). 

Proof of Theorem 14.81 Note that the second equality in (|4.26p follows from Theorem 
13.21 and Lemma 14.31 and similarly to the proof in Theorem 14.71 we only need to show that 
E{M,f ( J)M^{ J')} is super-exponentially small when J and J' are neighboring. Let I := 
J n J' 7^ 0. We follow the assumptions in (14.160 and assume also that all elements in e(J) 
and e(J') are 1, which implies E{J) = R^~ k and E(J') = R^~ k ' . 

We first consider the case k > 1. By the Kac-Rice metatheorem, E{M U ( J)M U (J')} is 
bounded from above by 



dt ds dx dy dz k+1 • • • / dz k+k ,_ m / dw m+1 ■ l dw k 
j J J 1 Ju Ju Jo Jo Jo Jo 

E{|detV 2 X| J (t)||detV 2 X| <7 ,( S )||X(t) = ar,X( S )=2/,VX| J (t) = 0,X fe+1 (t) = z fe+1 , 

. . .,X k +k'-m(t) = Zfc+fc'-m,VA|j/(s) = 0, X m+1 (s) = w m+ i, . . . ,X k (s) = w k } (4.27) 

X Pt,s( X > y> °> • • • ' z k+k'-m, 0, W m +l, ...,W k ) 

■ = A(t,s)dtds, 

J JjxJ' 
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where p t , s (x, y, 0, z k +i, . . . , z k +k'-m, 0, w m +i, ■ ■ ■ , w k ) is the density of 

(X(t), X(s), VX|j(t), X k+l (t), X k+kl _ m (t),VX l j / (s),X m+1 (s), X k (s)) 

evaluated at (x, y, 0, z k +i, • • • , Zk+k'-m, 0, w m+1 ,...,w k ). 

Let {ei, e2, • • • , eAr} be the standard orthonormal basis of M. N . For t 6 J and s G J', let 
et,s = (s - *) T /IN - t\\ and let a»(i,s) = (e i5 (A - A(i))e t)S ), then 

iV iV 

(A - A(t))et, a = J] (e^ (A - A(i))e M ) ei = ^ a<(t, (4.28) 

i=l i=l 
By Lemma 14.91 there exists some ao > such that 

(e M ,(A-A(i))e M ) > a (4.29) 

for all t and s. Under the assumptions (|4.16j) and that all elements in e(J) and e(J') are 1, 
we have the following representation, 

t = (ti, . . . , t m , i m +i, . . . , t k , b k +i, • • • , bk+k'-m, 0, . . . , 0), 
s = (si, . . . , s m , 6 m +i, . . . , b k , Sfc+i, . . . , Sfc + fc/_ m , 0, . . . , 0), 

where t, G (oj,6j) for all i G cr(J) and Sj G (a,j,bj) for all j G cr(J'). Therefore, 

(ei, e tyS ) > 0, V to + 1 < i < /c, 

(ei,e t)S )<0, V/c + l^i^/c + fc'-TO, (4.30) 
(ei,e M ) = 0, VHfc'-m<KJV. 



Let 



= {(t, s) G J x J' : ati(t, s) > Pi}, if m + 1 < i < k, 
Di = {(t, s) G J X J' : oti(t, s) < -Pi}, if k + 1 < % < k + k' - m, 

D = Ut,s) G J x J' : ^2ai(t,s)(ei,e t , s ) > Po>, 
^ i=l ' 



(4.31) 



where Po, Pi, ■ ■ - , Pk+k'-m are positive constants such that /3o + X^m+i" 1 ft < a o- It follows 
from (|4.30p and (|4.3ip that, if (i, s) does not belong to any of Do, D m , . . . , D k+k r_ m , then by 

(OH}, 

N k+k'-m 

((A - A(t))et, a , &t,s) = y^ajjt, s){ej, e t , s ) < Po+ ^ ft < a , 

i=l i=m+l 
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which contradicts P~29|) . Thus D U U^' +1 m A is a covering of J x J', by (f¥T27|) . 

E{Mf (J)Mf (J')} < / / A(t,s)dtds+ V / / A(t,s)dtds. 

J J D i=m+1 J J Dl 

We first show that j j Do A(t, s) dtds is super-exponentially small. Similarly to the proof 
of Theorem 113 applying (j4T20|) . P~2lj) and P~25j) . we obtain 

/ / s) cftfis 

< / / (iitis 



dxp V x, J (t) ) vx IJ ,( s )(0,0)p x(t) (a:|VX| J (*) = 0,VA'|j/(s) = 0) 
'Do A (4.32) 

x E{|detV 2 X| J (i)||detV 2 X| i/ /(s)||X(t) = x,VI[;(t) = 0,VX|j,(s) = 0} 
<C[ dtds dx(l + x N ')\\s-t\\- m p x{t) (x\VX lJ (t) = 0,VX|j/(s) = 0), 

for some positive constants C{ and iVi. Due to Lemma [4. 6\ we only need to consider the case 
when ||s — t\\ is small. It follows from Taylor's formula (|4.22p that as ||s — t\\ — > 0, 

Var(A-(t)|VX|j(t), VX {J ,{s)) < Var(A-(t)|^i(t), . . . , X m (i), X^s), . . .,X m (s)) 
= Var(X(t)|X!(t), . . . ,X m (t),Xx(t) + (VXi(t), s - i> + ||* - i|| 1+r? ^ s , . . . , 

X m (t) + (VJ m (t), s - t) + \\s - t\\ 1+r >Y t m s ) 
= Var(X(*)|*i(i), . . . ,X m (i), (VXi(t),e t)S ) + ||* - t\\W t ] s , . . . , (4.33) 

(VX m (t),e M ) + || S -t|TO) 
< Var(X(t)|(VXi(t), e t , s ) + ||s - ify£, . . . , (VX m (t), e t , s ) + ||s - ifY£) 
= Var(X(t)|(VX 1 (t), e M ), . . . , (VX m (t), e t)S » + o(l). 

By Lemma [6.2| the eigenvalues of [Cov((VXi(i), et )S ), . . . , (VX m (i), e^))] -1 are bounded 
uniformly in t and s. Note that E{X(i)(VXj(i), et jS )} = — aj(i, s). Applying these facts and 
Lemma 16.11 to the last line of (|4.33p . we see that there exist constants C2 > and £0 > 
such that for \\s — t\\ sufficiently small, 

m 

Y^{X{t)\VX\j{t),VX\j,{s)) <4-C 2 J> 2 (i, S ) + o(l) <4-£o, (4.34) 

i=i 

where the last inequality is due to the fact that (t, s) 6 Z?o implies 

Y> 2 (t,s) > Va 2 (t, S )|( ei ,e M )| 2 > - ( V on(t, s)(e t>s , ej) ) > ^. 
i=l i=l v i=l 7 
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Plugging (|4.34p into (|4.32p and noting that l/||s — t\\ m is integrable on J x J', we conclude 
that J J Do A(t, s) dtds is finite and super-exponentially small. 

Next we show that J f D , A(t, s) dtds is super-exponentially small for i = m + 1, . . . , k. It 
follows from (|4.27p that f J D A(t, s) dtds is bounded by 



oo roo 



dtds / dx / dwip x , t)t v X a )jX ., s}jVX , s) (x, 0,^,0) 

Di J a JO 

x E{|detV 2 X|j(t)||detV 2 X|j,(s)||X(t) = x,VX\j{t) = 0,Xi{s) = w^VX^^s) = 0}. 

(4.35) 

We can write 

Px(t),Xi(s){x,m\Xi{t) = 0) 



2ira 1 (t)a 2 (t,s)(l-p 2 (t 1 s)y/ 2 

{ 1 fx 2 wf 2p(t, s)xw,i 



where 



2(l-p 2 (i,s)) \a 2 {t) a 2 (t,s) a x {t)a 2 {t, a) ) J' 



aj(t) = V a r(X(t)\X l (t)=0), p(t,s) = - X ■' 1 ' 



<7l(i)cr 2 (t, s 

al(M) = Var(X i ( S )|X i (t) = 0) - ^(^(^^(t)) 



An 

and p 2 (t,s) < 1 due to (H3). Therefore, 

PX(t),VX| J (t),X i (s),VX| J /(s)( 2; i °' ^i) °) 

= Pvx } j l (s),x 1 (t),...,x i - 1 {t),x i+1 (t),,,,,x k {t)(S>\X(t) = x,Xi(s) = Wi,Xi(t) = 0) 

x Px{t),Xi(s)(x,Wi\Xi(t) = 0)p x . (t) (0) ( 4>36 ) 
f 1 f x 2 wf 2p(t, s)xwi \ ~) 

- 3 exp i 2(1 - P 2 (t, sj) v^!M + " *y j 

x (detCov(X(t), VX, VX, j,(s)))~ 1/2 
for some positive constant C3. Also, by similar argument in the proof of Theorem 14.71 there 
exist positive constants C4, C5, Cq, C7, ./V2 and iV"3 such that 

detCov(VX\j(t),Xi(8),VX\ jr (8)) > C A \\s - t\\ 2{m+1) , (4.37) 

C 5 \\s - t\\ 2 < a 2 2 (t, s) < C 6 ||s - t|| 2 , (4.38) 

and 

E{|detV 2 X|j(t)||detV 2 X|j,( S )||X(t) = x,VX\j(t) = 0,Xi(s) =w i ,VX\ JI (s) = 0} 
= E{\detV 2 X lJ (t)\\detV 2 X lJ ,(s)\\X(t)=x,VX l j(t) = 0, 

(4.39) 

(VIi(i),e t , s ) =^/||s-t|| + o(l),V*|j/(s) = 0} 
< 0-7(0^ + -tlD^ + l). 
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Combining (|4.35p with (|4.36p . (|4.37p and (|4.39p . and making change of variable w = Wi/\\s ■ 
t\\, we obtain that for some positive constant Cg, 



/ / A(t, s) dtds 

J JDi 



r r poo poo 

<C 8 dtds\\s - t|r m_1 / dx dwi(x N2 + (wi/\\s-t\\) N3 + 1) 

J JDi Ju JO 

j 1 fx 2 wf 2p(t, s)xwi \ 1 

x exp i - 2(1- P i(t,s)) [j*® + - ^rmJ / (440) 



P P POO PCO 

C 8 / / dtds\\s - t\\- m / dx dw(x N2 + w Na + 1) 
J Vd* Ju Jo 

( 1 / x 2 w 2 2p(t,s)xw 

xexpl - — -j— + 



where ^(t, s) = <J2{t,s)/\\s — t\\ is bounded by (|4.38j) . Applying Taylor's formula (14.221) to 
Xi(s) and noting that K{X(t)(VXi(t), et )S )} = — a^i, s), we obtain 

^' s ) = , * fE{X(t)X t ( S )} - -!-E{A:(t)X i (t)}E{X i ( a )X i (t)} 



s - t 



a 1 (t)a 2 (t,s) 



a i (t,s) + \\s-t\\ r >E{X(t)Y t ] s } 
\\s - t|K 



By (|4.38p and the fact that (t,s) G implies s) > > for i = m + l,...,k, 
we conclude that p(t,s) < —5q for some Jo > uniformly for t £ J, s £ J' with ||s — t\\ 
sufficiently small. Then applying Lemma f4. 101 to (|4.40p yields that f J D A(t, s) dtds is super- 
exponentially small. 

It is similar to prove that f f D _A(t,s)dtds is super-exponentially small for i = k + 
1, . . . , k + k' — m. In fact, in such case, J J D A(t, s) dtds is bounded by 



dtds 



oo coo 



dx I dz i px(t)yx l j(t),x l (t),vx lJ/ (s){x,0,z i ,0) 



x E{|detV^|j(t)||detV^|j*(a)||X(i) = x,VX\j(t) = 0,Xi(t) = z u VX VJ ,{s) = 0}. 

We can follow the proof in the previous stage by exchanging the positions of Xi(s) and Xi(t) 
and replacing Wi with Z{. The details are omitted since the procedure is very similar. 
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If k = 0, then m = and cr( J') = {1, . . . , k'}. Since J becomes a singleton, we may let 
J = {to}. By the Kac-Rice metatheorem, K{M U (J)M U (J')} is bounded by 



OO POO 



ds dx dy I dz 1 --- dz k , p to , s (x, y, z 1 , . . . , z k >, 0) 

Ju Ju JO JO 

xE{\detV 2 X lJ ,(s)\\X{t ) = x,X(s) = y,X 1 (t Q ) = z x , . . . ,X k ,{t Q ) = z h ,,V X\ r {s) = 0} 
A(t ,s) ds, 



where Pt , s (%, V, zi, . . . , z k >, 0) is the density of (X(i ), X(s), Xi(t ), • • • , -Xfc'(to), VX|j/(s)) 
evaluated at (x,y,z\, . . . ,z k ',0). Similarly, J' could be covered by U^ =1 Di with Dj = {s S 
J' : aj(to 5 s ) < — ft} for some positive constants ft, 1 < i < k'. On the other hand, 

r r roo /*oo 

_ A(t ,s)ds< ds dx dzip X ( to ^ Xt{to)yx {s) (x, Zi,0) 
Di Jd, Ju Jo 

x E{|detV 2 X|j,(s)||A:(io) = x,Xi{to) = z^VX^s) = 0}. 

By similar discussions, we obtain that E{M,f ( J)M^ ( J')} is super-exponentially small and 
hence complete the proof. □ 



5 Further Remarks and Examples 

Remark 5.1 [The case when T contains the origin] We now show that Theorem 14.71 and 
Theorem 14. 81 still hold when T contains the origin. In such case, (H3) is actually not satisfied 
since X(0) = is degenerate. However, we may construct a small open cube To contain- 
ing such that sup teTo of is sufficiently small, then according to the Borell-TIS inequality, 
P{sup tgTo X(t) > u} is super-exponentially small. Let T = T\Tq, then 

p{supX(i) > u) < p(supX(t) > u) < pjsupX(i) > u\ +P{ supX(t) > it). (5.1) 

To estimate P{sup te y > u}, similarly to the rectangle T, we decompose T into several 
faces by lower dimensions such that T = VJ^ =Q d k T = U^ L) LeQ ^ L. Then we can get the 
bounds similar to f|4. 3[) with T replaced with T and J replaced with L. Following the proof 
of Theorem 14.71 yields 

N 

P{ supX(i) > u} = Yl E { M «( L )i + o(e- au2 - u2/( - 2c7 T)). 
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Due to the fact that sup fgTo of is sufficiently small, K{M U (L)} are super-exponentially small 
for all faces L such that L C O^Tq with < k < N — 1 (note that To is a compact rectangle). 
The same reason yields that for 1 < k < N, L G d^T, J G such that L C J, the difference 
between E{M U (L)} and E{M U (J)} is super-exponentially small. Hence we obtain 

N 1 

P{supx (t )>„}= £ E , 

* 6T {t}69 r * fe=uea fe T v ^ 1 J| (5.2) 



Here, by convention, if Ot = 0, we regard e~ u2 ^ 29 ^ as 0. Combining (|5.ip with (|5.2|) . we 
conclude that Theorem 14. 71 still holds when T contains the origin. The argument for Theorem 
IE! is similar. 



Remark 5.2 Based on the proofs of Theorems 14.71 and 14.8^ one may expect that the 
approximation holds for a much wider class of smooth Gaussian fields (not necessarily 
with stationary increments). Meanwhile, the argument for the parameter set could go far 
beyond the rectangle case. These further developments will be included in Cheng (2013). 

Example 5.3 [Refinements of Theorem l4.7j Let X be a Gaussian random field as in Theorem 
14.71 Suppose that z^(to) = a\ for some to G J G d^T (k > 0) and v(t) < a\ for all t / to- 

(i). If k = 0, then, due to (|4. 14[) . sup tgT y| to } Q\ < a\ — Eq for some Eq > 0. This implies 
that E{M n (J')} is super-exponentially small for every face J' other than {to}- Therefore, 

P{ supX(t) >u} = *(—) + (e~ u2 /(24)+<*« 2 ), as u -)• oo. (5.3) 

I teT ) \(T T J 



For example, let Y be a stationary Gaussian field with covariance p(t) = e™"*" 2 and define 
X(t) = Y(t) — ^(0), then X is a smooth Gaussian field with stationary increments satisfying 
conditions (H1)-(H3). Let T = [0, 1]^, then we can apply (|5.3p to approximate the excursion 
probability of X with to = (1,...,1). 

(ii). If k > 1, then similarly, E{M U (J')} is super-exponentially small for every face 
J' 7^ J. It follows from Theorem 14.71 that 

Let r(t) = #f , then Vi G cr(J), Tj(to) = 0, since to is a local maximum point of r restricted 
on J. Assume additionally that the Hessian matrix 

ej(t ) := Mto))ij e „(j) (5-4) 
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(5.5) 



is negative definite, then the Hessian matrix of l/(2#f) at to restricted on J, 

©j(*0) = ~ 2T 2( t ^ ( T ij( t o))i,j£v(J) = ~2o T0J ^°' ) ' 

is positive definite. Let g(t) = \Aj - Aj(t) | /6<f fc ~ 1 and h(t) = 1/(26*^), applying Lemma [67 
with T replaced with J gives us that as u — > oo, 

ro f vu\^ \ «* -1 |Aj-Aj(t )| (2vr) fc / 2 . 

F< sup A m > It > = ; 7- — ^ kr-t e /V 'OMl + Oll)) 

Wer J (2vr)( fc +i)/ 2 |A J |i/2^-i M A:|e J ( io )|i/2 v v " 

2 fc / 2 |A J -A J fa)| , m . 

Example 15.31 [Continued: the cosine field] We consider the cosine random field on M 2 : 

1 2 

Z(t) = -=J2& COSt i + & s[nt i)> t=(tl,t 2 ) £ R 2 , 
v 2 i= i 

where £i, £ 2 , £2 are independent, standard Gaussian variables. Z is a well-known centered, 
unit- variance and smooth stationary Gaussian field [cf. Adler and Taylor (2007, p. 382)]. 
Note that Z is periodic and Z(t) = —Zu(t) — Z 22 (i). To avoid such degeneracy, let X(t) = 
£0 + Z(t) - Z(0), where t £ T C [0, 2vr) 2 and £ is a standard Gaussian variable independent 
of Z. Then X is a centered and smooth Gaussian field with stationary increments. The 
variance and covariance of X are given respectively by 

u(t) = a t = 3 — costi — cos £2, 

(5.6) 



C(t,s) = 2 H — y^[cos(tj — Si) — costj 



, COS.S,|. 

8=1 



Therefore, A" satisfies conditions (HI), (H2) and (H3) [though Xi 2 {t) = 0, it can be shown 
that this does not affect the validity of Theorems 14.71 and I4.8j . Taking the partial derivatives 
of C gives us that 



E{X(t)VX(t)} = ^(sinii,sint 2 ) T , A = Cov(VX(t)) = h 2 , 

A — A(i) = -E{AT(t)V 2 X(i)} = -[J 2 -diag(costi,cosi 2 )], 



(5.7) 



where I 2 is the 2x2 unit matrix and diag denotes the diagonal matrix. 

(i). Let T = [0,7r/2] 2 . Then by (|5.6h . v attains its maximum 3 only at the corner 
(tt/2, 7r/2), where both partial derivatives of v are positive. Applying the result (i) in Example 
IQl we obtain P{sup teT X(t) > u} = (1 + o(e~ att2 )). 



28 



(ii) . Let T = [0, 37r/2] x [0,7r/2]. Then v attains its maximum 4 only at the boundary 
point t* = (7r,7r/2), where ^ 2 (i*) > so that the condition (|4.14|) is satisfied. In this case, 
t* G J = (0,3vr/2) x {vr/2}. By ([573), we obtain Aj = \ and Aj-Aj(t*) = ±(1- cos if) = 1. 
On the other hand, for t G J, by Lemma UTTl and (|5.7|) . 

r(t) = 6» 2 = Var(X(t)|Xi(t)) = 3 - cosii - cosi 2 - ^sin 2 £i, (5.8) 

therefore Oj(i*) = Tu(t*) = —2. Plugging these into (15. 5D with = 1 gives us that 
P{sup ieT X(t) >u} = >/2¥(u/2)(l + o(l)). 

(iii) . Let T = [0, 37r/2] 2 . Then i/ attains its maximum 5 only at the interior point t* = 
(vr,7r). In this case, t* £ J = (0,3tt/2) 2 . By ifBTFj) . we obtain Aj = |l 2 and Aj - Aj(i*) = I 2 . 
On the other hand, for t G J, by Lemma 16. II and (|5.7p . 

r(i) = (9 2 = Vax(X(t)\Xi(t),X 2 (t)) = 3 - costi - cost 2 - ^sin 2 ti - ^sin 2 t 2 , (5.9) 

therefore Qj(t*) = (t^- j=i,2 = — 2/ 2 . Plugging these into (15. 5p with k = 2 gives us that 
P{sup teT X(t) >u} = 2^(u/ v / 5)(l + o(l)). 



Example 5.4 [Refinements of Theorem 14. 8j Let X be a Gaussian field as in Theorem 14.81 
Suppose to £ J £ dkT is the only point in T such that v{to) = o\. Assume cx( J) = {1, . . . , k}, 
all elements in e(J) are 1, Vk'(to) = for all k + 1 < k' < N. Then by Theorem 14.81 

N 

P{ sup X(t) > = E{Aff (J)} + E{Mf (J')} + o(e- au2 -" 2/(2<T T)). 

" fc'=fc+i j'GO fc /T,J'nJ^0 

(5.10) 

Lemma HTJ indicates E{Mf (J)} = (-l) fc E{ELo(- 1 )Vi(^)}(l + o(e" aa;2 )), therefore 



E{Mf(J)} = (-l) fc / pv^^^WdiEldetV 2 ^^*)^^ 



(i),...,x JV (t))eR^-' = } 



x V(*)>«}I VA V(*) = 0}(1 + oieT** )) 



{(X fe+1 (t),...,X JV (t))GR^- fc ] 



|X(t) = x,VA| J (t) = 0}(l + O (e- a « 2 )) 
^j(x)dx(l + o( e - Q " 2 )), 



(5.11) 



and similarly, 

, (_Dfc' e -^ 2 /(2e t 2 ) 



E{M-(J0} = / y dt n^ eW^At) 



j, (2vr)( fc '+ 1 )/ 2 |A J /| 1 /2^ 

x \ Xjl{t) _ Xjl m& »-« } W) = *,vx\At) = o}(i + o( c — 2 )). 
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(i). First we consider the case k > 1. We shall follow the notation r(i), Qj(t) and Qj(t) 
in Example E3 Let h(t) = 1/(26*^) and 

( — l) fc 2 

= (2vr)(^i)/2|A J |i/20 t E{detV2X l J ^ )1 {(^ + iWv..,x i v( t )) G Rrn 

\X(t) = x,VX ]J (t) = 0}. 

Note that sup tgT Ifl^)! = ^(x^ 1 ) for some N% > as x — > oo, which implies that the growth 
of g x (t) can be dominated by the exponential decay e~ x2h ^\ hence both Lemma 16.31 and 16.41 
are still applicable. Applying Lemma 16.31 with T replaced by J and u replaced by x 2 , we 
obtain that as x — > oo, 

Aj(x) = —^§^77 1/2 fe(to)e- 2 ^(l + (5.12) 
x fe (detej(i )) i/i! 

On the other hand, it follows from (|3.17p that 

*>® = w + iyz\Aj\v*e t I • • • / Kr , dVk+l - dm 

x |Aj ~f j(t)l g fc ( - + 7tCfc+i(t)yfe+i + • • • + ltC N (t)y N 
It V7t 

x Px fc+1 (t),...,x w (t)(yfc+i, • • .,2/jv|X(t) = x,VX\j(t) = 0). 

Note that X(to) and VX(to) are independent, and Cj(to) = for all 1 < j < N. Therefore, 

s _ \Aj-Aj{t )\ fx 
9x[ 0) ~ (2^)(^)/2|A J |i/24+ 1 

x P{(X fc+ i(to), . . . ,X N (t Q )) G Rj-*|VX|j(i ) = 0}. 
Plugging this and (I5.12p into (j5.11j) . we obtain 

E{M E (J)} = 2k/2 ^~ A ^\ 

x ¥{(X k+1 (t ), . . .,X N (t )) G Ry- fc |VX|j(* ) = 0}(1 + o(l)). 

For J' G <9fc/T with J' n J 7^ 0, similarly, applying Lemma 16.41 with T replaced by J', we 
obtain that 

«wf (-0) = rfw ^ <-)n^(^o) € 

|Aj'| v I - ©J'^o)! 172 v o- T / (5.14) 
x P{(X j; (t ), • • • ,Xr .(to)) e R^'lVX^ito) = 0}(1 + o(l)), 

1 N — k' 1 
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where Zjt(ta) is a centered (k' — /c)-dimensional Gaussian vector with covariance matrix 
~( T ij)i,jecr(J')\<T{J)- Plugging (|5.13|) and (|5.14p into (|5.1U|) . we obtain the asymptotic result, 
(ii). k = 0, say J = {to}- Note that X(t ) and VX(i ) are independent, therefore 

E{Mf (J)} = (-^-)p{VX(t ) E M+}. (5.15) 

For J' G <9 fc ,T with J'n J ^ 0, then E{M^(J')} is given by ([5T4]) with /c = 0. Plugging ([5T5]) 
and (|5.14p into (|5.10p . we obtain the asymptotic formula for the excursion probability. 
Example 15.41 [Continued: the cosine field] We consider the Gaussian field X defined in 
the continued part of Example 15.31 

(i) . Let T = [0,7r] 2 . Then v attains its maximum 5 only at the corner t* = (vr, 7r), where 
Vu{t*) = so that the condition (|4.14p is not satisfied. Instead, we will use the result (ii) in 
Example with J = {t*} and k = 0. Let J' = (0, vr) x {vr}, J" = {vr} x (0,vr). Combining 
the results in the continued part of Example 15.31 with (|5.15p and (|5.14|) , and noting that 
A = J? -£2 implies X±(t) and X%(t) are independent for all t, we obtain 

E{Mf (J)} = ^(u/x/5), E{Mf (fcT)} = ^*(u/VS)(l + o(l)), 
E{Mf (J')} = E{Mf (J")} = ^*(u/V5)(l + o(l)). 

Summing these up, we have P{sup teT X(t) > u} = [(3 + 2 v / 2)/4]*(n/ v / 5)(l + o(l)). 

(ii) . Let T = [0, 37r/2] x [0, ir]. Then v attains its maximum 5 only at the boundary point 
t* = (-7T, 7r), where v 2 {t*) = 0. Applying the result (i) in Example with J = (0, 37r/2) x {7t} 
and k = 1, we obtain 

E{Mf (J)} = ^¥(«/V5), E{Mf (0 2 T)} = *(ti/V5)(l + o(l)), 
which implies P{sup tgT X(t) > u} = [(2 + \/2)/2]*(u/V5)(l + o(l)). 

Remark 5.5 Note that we have only provided the first-order approximations for the exam- 
ples in this section. However, as shown in the theory of approximations of integrals [see e.g. 
Wong (2001)], the integrals in (|4.15p and (|4.26p can be expanded with more terms provided 
the covariance function of the Gaussian field is smooth enough. Hence for the examples above, 
higher-order approximation is available. Since the procedure is similar and the computation 
is tedious, we omit such argument here. 

6 Appendix 

This Appendix contains proofs of Lemmas 14.11 14.91 14.101 and some other auxiliary facts. 
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Proof of Lemma 14.11 By the definition of (J), it is clear that 

TV 

X(t) > u] D 



iV 

{ supX(t) >ii}d[J (J {Mf (J) > 1} a.s. 



fc=o Jed fc T 

Suppose sup tgT A"(i) > u, since X(t) G C^IR^) a.s., there exists to £ T such that X(fo) = 
sup tgT X(i). Without loss of generality, assume to £ J £ dkT. Note that to is a local 
maximum restricted on J, thus VX\j{to) = and V 2 X\j(to) is non-positive definite. Due 
to (HI) and (H3'), we apply Lemma 11.2.11 in Adler and Taylor (2007) to obtain that 
almost surely, det(V 2 X| j(t )) ^ and hence index(V 2 X| j(t )) = k. U e*Xj(t ) < for some 
j ^ &(J), then we can find t\ £ T such that X(ti) > X(to), which contradicts X(to) = 
sup teT X(t). Hence e*Xj(t ) > for all j g a (J). These indicate (J) > 1, therefore 



N 

\ S upX(t) > v) C M I) {Mf (J) > 1} a.s., 



fc=0 Je9 fc T 

completing the proof. □ 

Proof of Lemma 14.91 Let M^vxTV be the set of all N x N matrices. Define a mapping 
(f) : M N x MnxN — > K by (£,^4) h-> (£, A£), then is continuous. Since A — A(i) is positive 
definite, 0(e, A - A(t)) > for each i G T and e G S^ -1 . On the other hand, {(e, A - A(i)) : 
t £ T,e £ S^ -1 } is a compact subset of K N x M^vx/v and is continuous, completing the 
proof. □ 

Proof of Lemma 14.101 We only prove the first case, since the second case follows from the 
first one. By elementary computation on the joint density of and ^(s), we obtain 

sup E{(1 + |&(t)|^ + \Us)\ N2 n { ^t)>uMs)<o}} 



2 



u 2 gjp 2 ^ 2 2 

o exp < ^ — ^ „, - — — + eu 

\ l \ 2a\ 2a 2 (l-p 2 )a 2 

as u — > oo, for any e > 0. □ 
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The following lemma is well-known and is quoted here for reader's convenience. 



Lemma 6.1 LetY and Z be two Gaussian random vectors of dimension p andq, respectively. 
Then Y\Z = z is still a p-dimensional Gaussian random vector having the following mean 
and covariance: 

E{Y\Z = z} = EY + E{(Y -EY)(Z - EZ) T }[Cov(Z)] _1 (z -EZ), 
Cov(Y\Z = z) = Cov(Y) - E{(Y - EY)(Z - EZ) T }[Cov(Z)] _1 E{(Z - EZ)(Y - EY) T }. 

In particular, if p = q = 1 and EY = EZ = 0, then 

Using similar argument in the proof of Lemma |4.9| we can obtain the following result. 

Lemma 6.2 Let {A(t) = (oy (t))i<ij<iv ■ t £ T} be a family of positive definite matrices 
such that all elements are continuous. Denote by x and x the infimum and supremum 

of the eigenvalues of A(t) over t £ T respectively, then < x < x < oo. 

The following two formulas state the results on the Laplace method approximation. 
Lemma 16.31 can be found in many books on the approximations of integrals, here we re- 
fer to Wong (2001). Lemma 16.41 can be derived by following similar argument in the proof of 
Laplace method for the case of boundary point in Wong (2001) [see Cheng (2013)]. 

Lemma 6.3 [Laplace method for interior point] Let to be an interior point of T. Suppose 
the following conditions hold: (i) g(t) G C(T) and g(t$) ^ 0; (ii) h(t) G C 2 (T) and attains 
its unique minimum at t$; and (Hi) V 2 /i(fo) is positive definite. Then as u — >■ oo, 

//W'^''- tt = ^ffl))^ « ( '° ) '^ (W(1+0(1)) - 

Lemma 6.4 [Laplace method for boundary point] Let to G J G d^T with < k < N — 1. 
Suppose that conditions (i), (ii) and (Hi) in Lemma \6.S\ hold, and additionally V/i(io) = 0. 
Then as u -4 oo, 

where Zj(to) is a centered (N — k)- dimensional Gaussian vector with covariance matrix 
(hij(to))j 1 <ij<j N _ k , —E(J) = {x G : — x G E(J)}, and the definitions of J\, . . . , Jjv-fc 
and E(J) are in \3.^- 
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